£(1385) Results and Status of the 9+ in STAR 



Sevil Salurf for STAR Collaboration! 

Physics Department, Yale University, 
Sloane Physics Laboratory, P.O. Box 208120, 
New Haven, CT, 06520-8120, USA 



O 



(j • Abstract. The £(1385) analysis and the current status of pentaquark search with 

, the STAR detector are reported. The corrected px spectra and the yields of the 

E ± (1385) and their antiparticles in the most central Au+Au as well as elementary 
• P+P collisions are presented. A comparison of the (px) of observed particles suggests 

a similar behavior for particles with mass greater than 1.2 GeV in p+p and Au+Au 
collision environments. Acceptance and efficiency studies with simulations show that 
. the (anti)pentaquarks should be found at the 3 % level. 
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■^j- ■ 1. Introduction 

o ' 

The observation of a five-quark bound system consisting of uudds, referred to as the B + 
pentaquark, have been reported in photon-nucleus and kaon-nucleus reactions [1, 2, 3]. 
The presence of this state was predicted by R. L. Jaffe with multiquark bag models 
[4, 5] and later by D. Diakonov et al. using chiral soliton models of baryons [6]. The 
high energies and particle densities resulting from collisions at the Relativistic Heavy 
Ion Collider (RHIC) are expected to be ideal environments for pentaquark production 
[7, 8, 9, 10]. 

During the expansion of the hot and dense matter (fireball) created in Au+Au 
collisions, chemical freeze-out is reached when the hadrons stop interacting inelastically. 
Elastic interactions continue until thermal freeze-out. Due to the very short lifetime 
(r < r fireball) of resonances, a large fraction of the decays occur during this time. 
The elastic interactions of decay products with the surrounding particles result in a 
signal loss in the particle identification though this is offset by secondary interactions 
which increase the resonance yield (e.g. A + 7r — ► £(1385), K+p — > @ + ). The 
contribution of re-scattering and regeneration on the total observed yields depends on 
the time span between the chemical and thermal freeze-out and the lifetime of each 
resonances [11, 12]. Thus the study of resonances provides an additional tool in the 
determination of the hadronic expansion time between chemical and thermal freeze-out 
by comparing resonance to stable particle ratios. The data analyzed were taken by the 
STAR (Solenoidal Tracker At RHIC) experiment [13], one of the four experiments at 
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RHIC. The large acceptance of STAR's Time Projection Chamber (TPC) is ideal for 
such rare particle searches. 

2. Particle Identification 

The STAR trigger detectors consist of two zero degree calorimeters (ZDCs) which are 
situated ~ 17 m downstream of the nominal interaction point as well as a central 
trigger barrel (CTB) surrounding the TPC. A minimum bias trigger is defined by the 
coincident measurement of spectator neutrons in both ZDCs. The CTB an array of 
scintillator slats, is used to detect event multiplicity is used to trigger on the 10% most 
central Au+Au collisions. For p+p collisions, a minimum bias trigger is defined by 
coincidences in the two beam-beam counters, which are again scintillator detectors and 
are situated around the beam pipe, approximately 2 m from the center of the TPC. The 
main detector component of STAR is the TPC, which, together with the magnetic field 
information, is used to identify stable charged particles via energy loss per unit length 
and long-lived weakly decaying (cr ~ few cm) neutral particles such as A and K s via 
their decay topology. The direct measurement of resonances is not possible due to their 
short lifetimes (cte(i385) = 5 fm) [14, 15]. Instead, £(1385) resonances are identified 
by the invariant mass after combining n with A decay particle candidates. A mixed 
event technique, where particles from different events are reconstructed using the same 
technique, is used to determine the background for uncorrelated pair combinations. The 
£(1385) signal is obtained by the subtraction of this normalized mixed event background 
from the invariant mass distribution [16] . This technique can directly be applied for the 
search and possible identification of pentaquarks. For the case of the 6 + pentaquark 
study the invariant mass reconstruction is via 6 + — > K s + p, and subtraction of the 
normalized mixed event background is used [18]. 

3. The £(1385) Analysis 

The transverse mass (% = \/pj> + m2 ) spectra of £ ± (1385) and their antiparticles 
in p+p (circles) and Au+Au (stars) collisions are shown in Fig. 1. The spectra are 
corrected for acceptance and efficiency by embedding Monte-Carlo simulated resonances 
into real p+p and Au+Au events. The solid lines in Fig. 1 represent exponential fits 
to the data with the function directly proportional to the yield (dN/dy) and inversely 
proportional to temperature (T 2 + m T). The data coverage is 91% in p+p and 85% 
in Au+Au collisions of the fully integrated yield. The (px) is derived from the full 
range integration of the corresponding exponential fit. Table 1 presents the inverse 
slope parameter (T), the (jpr) and the yield (dN/dy) for the summed signal of £ ± (1385) 
together with their antiparticles, for both p+p and Au+Au collisions at y/s NN = 200 
GeV. 

The (pt) for various particles in both p+p and Au+Au collisions at ^s^n = 
200 GeV as a function of their mass are presented in Fig. 2. The behavior of (p T ) 
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Figure 1. The transverse mass spectrum for £(1385) drawn as stars for the — 5% 
most central Au+Au and circles for p+p collisions at ^Snn = 200 GeV. Blue circles 
are for p+p and red stars are for Au+Au collisions. 



vs. mass for the various particles in p+p and Au+Au collisions is compared to two 
parameterizations. The black curve is an empirical fit to the ISR tt, K and p data [19] 
and the shaded band is a blast wave fit using tt, K and p data in Au+Au collision system 
[17]. The empirical parametrization for the ISR data at y/s = 25 GeV in p+p collisions, 
can describe the behavior of the lower mass particles, such as tt, K and p, despite the 
fact that our collision energy is one order of magnitude higher. However, this empirical 
parametrization does not represent the behavior of the higher mass particles. Similarly 
the blast wave parametrization which can describe the lower mass particles in Au+Au 
collisions (~ 98% of all the particles observed) fails to explain the behavior of higher 
mass particles. 

The comparison of lighter, weakly decaying particles and heavier resonances in 
p+p and Au+Au collision environments shows a similar behavior of (pt) for the higher 





T [MeV] 


(Pt) [GeV/c] 


Yields (dN/dy) 


£±(1385) in p+p 


358 ± 47 


1.08 ±0.15 


(4.66 ±0.98) x 10~ 3 


£±(1385) in Au+Au 


420 ± 84 


1.20 ±0.24 


4.72 ± 1.38 



Table 1. Temperature T, (pt) and yield obtained from the exponential fits of 
the pt spectra in Fig. 1 for elementary p+p and 0-5% most central collisions. 
The statistical uncertainties are given and the systematical error due to norm- 
alization of the background, ~ 15%, has to be included in the given values. 



£(1385) Results and Status of the + in STAR 



4 



„ 1.6 r 
o 

> 
O 1.4 



1.2 
1 

0.8 

0.6- 

0.4 



s NN = 200 GeV 



s NN = 200 GeV 



% Central Au+Au 
MinBias p+p 

Blast Wave Parametrization 
=89+1 MeV, < P T > =0.59 ± 0.05 
ISR, < pr> = 0.7 X m' 

STAR Preliminary 



0.2 




1.4 1.6 1.8 
Mass [GeV/c 2 ] 



Figure 2. The (pr) vs particle mass measured in p+p and Au+Au collisions at 
\J snn — 200 GeV. While the circles represent the stable particles, triangles represent 
the resonances. The black curve represents the ISR parametrization from 7r, K and p 
for v'snn — 25 GeV p+p collisions. The shaded band is the blast wave fit using tt, K 
and p for Au+Au collisions. 



mass particles. It is expected that resonances with higher transverse momentum are 
more likely to be reconstructed in Au+Au collisions because of their longer relative 
lifetimes due to Lorentz contraction. This means they are more likely to decay outside 
the medium and hence their daughter particles would interact less with the medium. 
Any loss at low px would increase the T parameter of the px spectra for the central 
Au+Au collisions with respect to p+p collisions. However we do not see any significant 
increase in the T parameter for £(1385) from p+p to the most central Au+Au collisions 
within the statistical and systematic errors. If higher mass particles are produced by 
different production mechanisms than the lower mass particles, we might be introducing 
a bias in our measurement in p+p collisions. This may be due to a number of reasons 
which include higher mass particles being produced in more violent (mini-jet) p+p 
collisions and/or the heavier particles flow radially with a smaller velocity than the 
lighter mass particles (such as n mesons) in the most central Au+Au collisions. The 
(Pt) measurement of the higher mass resonance £(1385) shows this merging behavior. 
Similarly the study of radial flow of heavier particles can give us valuable information 
about whether the heavier particles flow less with respect to lower mass particles. 

A comparison of £(1385)/A ratios in p+p and — 5% central Au+Au collisions 
shows no suppression from p+p to Au+Au. An observed suppression of K*(892)/K 
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Figure 3. Invariant mass spectrum of the Monte Carlo simulated + embedded in 
real p+p events with the mixing technique. Black solid histogram is the signal and red 
dashed histogram is the mixed event background. The simulated signal can be clearly 
seen for 1 simulated particle per event. 



and A*(1520)/A ratios in Au+Au collisions with respect to p+p collisions suggests a 
signal loss via re-scattering in the medium by thermal models [20, 21]. As the lifetime 
of a £(1385) is one third that of the A*(1520), it is expected that there should be an 
even greater signal loss via re-scattering. The observed unsuppressed ratio requires a 
significant regeneration mechanism in order to recover the signal loss via re-scattering. 
The measured £ ± (1385)/A ratio is 0.295 ± 0.086 for the 0-5% most central collisions 
which is about a factor of 2 below the microscopic model UrQMD predictions at 
v /J/vtv = 200 GeV. This suggests that the assumed regeneration cross-section, included 
in UrQMD calculations, is too high [22]. Both UrQMD and thermal production models 
should be revised in light of the resonance measurements. 

4. Status of Current + Studies 

4-1. Monte Carlo Studies 

To study the decay mechanism and optimize the applied cuts, Monte Carlo simulations 
are used. In this study, one Monte Carlo G + pentaquark is chosen from a thermal 
exponential distribution with T = 250 MeV in the rapidity interval | y |< 1.5 and after 
a full TPC simulation is embedded into a single, real p+p event. The chosen input width, 
10 MeV/c 2 , and the input mass, 1.54 GeV/c 2 , are consistent with the observed mass 
and width of + [1, 2, 3]. In Fig. 3 the invariant mass spectrum of the reconstructed @ + 
is presented. We find that ~ 3% of these Monte Carlo generated 6 + 's are successfully 
reconstructed with this technique. The reconstructed width and the mass is consistent 
with the Monte Carlo input. 

Using this technique, the decay properties with the simulated tracks such as the 
momentum distribution of the decay daughters can be studied. In Fig. 4, the momentum 
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Figure 4. On the left is the Kg momentum distribution and on the right is the proton 
momentum distribution. Black solid histograms is of the accepted K s on the left and 
p on the right after the dE/dx cut and red dashed histograms is of the decay daughters 
of the Monte Carlo generated + for the same number of events. 



distributions of the Kg on the left and proton on the right are presented. With this study 
we can optimize momentum cuts to increase the signal-to-background ratio. Clearly, 
an optimized cut to improve this ratio is to accept protons with momentum less than 
1 GeV/c while at the same time rejecting the ones above this threshold. Further detailed 
studies on other variables are needed to optimize the signal-to-background ratio. 



4-2. Feasibility Studies 

Assuming that the B + production is 10 — 100% of the A* (1520) in p+p collisions, one 
can estimate the yield of the + . The preliminary dN/dy of A*(1520) at mid-rapidity 
is 0.004 per event in p+p collisions [23, 24]. There are 8 Million p+p events available 
for this analysis and this corresponds to a production of ~ 30000 A* (1520), giving a 
production range of ~ 3000 to ~ 30000 6 +, s in these p+p events using the above 
estimates. As the efficiency of the mixing technique is ~ 3% and the branching ratio of 
the + — > Kg + p is ~ 25% (assuming that the branching ratios of — > K s + p and 
+ — > K + +N are each 50%), 20-200 of the + 's should be found. The background pairs 
per event in the 1.54+5 MeV mass range is 3200. This corresponds to a significance of 
between 0.25 and 3 Similarly one can repeat the same study for Au+Au and d+Au 
collisions and correspondingly predict a significance of 2-7 for 1.5 Million Au+Au events 
and 1-16 for 10 Million d+Au for the predicted production of one + per unit rapidity 
per collision [7, 8, 9, 10]. To estimate the yield for the d+Au collisions we assume scaling 
with the number of participants (N part ). The mean number of participants in d+Au is 
8, in p+p it is 2, and in Au+Au it is 350 for the most central collisions. The lower limit 
is obtained from p+p scaling while the upper limit is from Au+Au yield estimates. The 
invariant mass spectra that are observed for the + in p+p, d+Au and Au+Au collision 
events are consistent with our estimations for the significance of the signal given our 

i The significance is defined as ,„ D Sl 9 nal 
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current statistics. The signal-to-background ratio depends highly on the selection of 
events and applied cuts. To improve cuts and understand the decay mechanism more 
detailed simulation studies must be undertaken. 

5. Conclusions 

£ ± (1385) resonances and their antiparticles can be reconstructed via event mixing 
techniques in p+p, d+Au and Au+Au collision environments with STAR. The £(1385) 
(pr) in p+p collisions is similar within the errors to that measured in Au+Au collisions. 
This follows the behavior of the other high mass particles (to > 1.2 GeV) whose (p?) 
measurement in p+p approaches the Au+Au value. This behavior suggests a possible 
smaller radial flow for heavy particles with respect to 7r and/or a more violent production 
mechanism for heavy particles in p+p collisions. 

There is no significant suppression observed in the £(1385)/A ratios from p+p 
to — 5% central Au+Au collisions. Since the previously measured suppression 
of K*(892)/K and A* (1520) /A ratios agrees with a signal loss via rescattering in 
the medium, the absence of suppression in the £(1385) ratio suggests a significant 
regeneration mechanism to recover the signal lost via rescattering [25, 26]. 

Acceptance and efficiency studies show that we should be able to find @ + 
pentaquarks at the few % level with the current data-set in d+Au and Au+Au collisions. 
Optimization of cuts to improve the signal over background is in progress. There is a 
possibility of measuring anti-pentaquarks at RHIC since the antibaryon to baryon ratio 
is approaching to one [27]. An upper limit to the yields and production mechanisms 
of pentaquarks in Au+Au collisions will be established in the 2004 run, which has a 
70-fold increase in statistics over the current data. 
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